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1.  Introduction 


The  combustion  synthesis  glass  using  chlorinated  silicon  compounds  (SiCl4/HSiCl3)  is  important  to 
industry.  The  AEROSIL®  process  is  a  well-established  process  for  generating  silica  from  high- 
temperature  combustion  of  chlorosilanes  in  flames.  Although  the  process  has  been  in  use  for  60 
years,  there  is  still  very  little  understanding  of  the  basic  underlying  chemical  reactions  involved  in  the 
process.  Currently,  very  little  is  known  about  the  SiCl4  oxidation  mechanism,  thermochemistry,  or 
kinetics  involved.  Several  studies  of  the  pyrolysis  of  chlorinated  silanes  have  recently  appeared.1-3 
Thermal  decomposition  of  SiCl4  has  been  measured  behind  reflected  shock  waves.2,3  SiCl3  was  found 
to  form  by  direct  thermal  decomposition  of  SiCl4.  ’ 


SiCl4 - >SiCl3  +C1 


(1) 


SiCl3  can  also  be  formed  by  the  reaction  of  SiCl4  with  H  in  the  gas  phase  at  high  temperatures.3,4 


H  +  SiCl4 - >SiCl3  +  HCl  (2) 

Recently,  reaction  (2)  was  found  to  be  a  major  initiation  reaction  in  SiCl4/H2/02/Ar  flames.4  Previous 
models  of  the  oxidation  process  have  assumed  spontaneous  formation  of  Si02  that  coagulate  to  form 
particles.5,6  A  model  that  includes  detailed  chemical  kinetics  has  recently  been  developed.4  The 
multi-step  reaction  mechanism  was  developed  by  modeling  experimental  observations  of  OH  forma¬ 
tion  in  a  SiCl4/H2/02  flame.4 

TiCl4  is  used  in  the  production  of  titania  and  TiN.7'8  TiCl4  can  be  added  to  SiCl4/H2/02  flames  to  pro¬ 
duce  Ti02  dopant  in  Si02  glass  and  ceramics.  Thermal  decomposition  of  TiCl4  has  been  measured 
behind  reflected  shock  waves.9  TiCl3  was  found  to  form  either  by  direct  thermal  decomposition  of 
TiCl4,  and  TiCl2  was  found  to  form  by  direct  thermal  decomposition  of  TiCl3. 


TiCl4  +  M- 

— *  TiCl3  +  Cl  +  M 

(3) 

TiCl3  +  M- 

— >  TiCl2  +  Cl  +  M 

(4) 

The  initiation  reaction  of  TiCl4  with  H  has  also  been  measured  behind  reflected  shock  waves  at  high 
temperature.10  There  is  little  kinetic  information  on  the  oxidation  of  these  TiClx  intermediates. 
Recently,  Fourier  transform  infrared  spectroscopy  was  used  to  observe  Ti02  formation  in  premixed 
flames  of  titanium  tetraisopropoxide,  methane,  and  oxygen.11  Previous  modeling  studies  have 
focused  on  the  coagulation  of  the  Ti02  particles.12,13  Teyssandier  and  Allendorf  have  previously 
modeled  the  Ti-Cl-H  system  by  using  Rice-Ramsberger-Kassel-Marcus  to  predict  reaction  rates.14  In 
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order  to  do  more  detailed  modeling  of  TiCU  oxidation  in  flames,  the  basic  kinetic  and  thermodynamic 
information  is  needed  on  intermediate  Ti-O-Cl-H  species. 

Theoretical  methods  have  previously  been  used  to  calculate  chlorinated  silane  thermochemistry.15'23 
Melius  et  al.  have  recently  used  M0ller-Plesset  perturbation  theory  with  bond  additivity  corrections 
(BAC-MP2  or  BAC-MP4)  to  predict  the  thermochemistry  of  several  silicon-containing  systems  (Si- 
C-O-H,  Si-O-H,  Si-Cl-H,  Si-H-O-Cl).15’16’21"23  The  BAC-MP4  method  is  believe  to  have  an  accuracy 
of  (±  2-7  kcal  mol-1).21  The  bond  additivity  correction  (BAC)  method  uses  an  empirical  method  to 
add  “corrections”  to  the  ab  initio  calculated  heats  of  formation.  In  order  to  get  good  empirical  cor¬ 
rections,  a  good  subgroup  of  experimentally  measured  species  would  be  needed  for  comparison. 
However,  for  the  titanium  flame  intermediate  species  of  interest,  very  little  thermochemistry  is  cur¬ 
rently  known. 

This  work  uses  the  MP2  ab  initio  method  to  calculate  several  postulated  intermediate  species  associ¬ 
ated  with  TiCU/SiCU/TB/Ch  flames.  The  error  associated  with  using  the  MP2  method  is  compared  to 
the  error  in  previous  theoretical  calculations  and  experiments. 


2 


2.  Theoretical  Method 


Electronic  structure,  normal  mode  analysis,  and  heats  of  formation  were  calculated  by  using  the 
Gaussian  03W24  ab  initio  program.  The  electronic  structure  was  first  optimized  by  using  the  unre¬ 
stricted  Hartree-Focke  (HF)  method  with  6-31G(d)  basis  set.  The  optimized  structure  was  then  used 
to  calculate  (HF/6-31G(d))  the  vibrational  frequencies  and  thermal  correction  to  enthalpy  (which  was 
scaled  by  0.8929).  The  electronic  structure  was  then  refined  at  the  MP2  level  of  theory  using  the  6- 
31G(d)  basis  set.  The  electronic  energy  was  then  calculated  using  MP2  with  6-31G(3df,2p)  basis  set. 
The  AHf°  (298K)  were  calculated  from  the  electronic  energy  (E)  of  the  species  of  interest  by  Eqs.  3-5. 

AE  =  E(SiaClbOcHdTie )  -  a  *  E(Si)  -  b  *  E(C1)  -  c  *  E(O)  -  d  *  E(H)  -  e  *  E(Ti)  (3) 

H  =  AE  +  Et (SiaClbOcHdTie )  -  a  *  ET(Si)  -  b  *  ET (Cl)  -  c  *  Ex (O)  -  d  *  ET (H)  -  e  *  ET  (Ti)  (4) 

AH° (0  K)  =  H-a*  AH°  (Sia ) - b *  AH° (Cl) -c*  AH°(0)-d  *  AH°  (Hd)-e*  AH° (Ti)  (5) 

where  a,  b,  c,  e,  and  d  are  coefficients,  and  ET  is  the  calculated  thermal  correction  to  the  electronic 
energy.  The  heats  of  formation  for  the  elements  where  taken  from  the  NIST  database  and  are  listed  in 
Table  l.25 


Table  1 .  The  Heat  of  Formation  of  the  Elements25  Used  to 
Calculate  the  Heats  of  Formation  Using  Eqs.  3-5 


Element 

AH°  (298K)  kJ  mol"1 

Si 

450. 

Cl 

121.30 

H 

218.00 

0 

249.17 

Ti 

473.63 
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3.  Results  and  Discussion 


The  calculated  MP2  electronic  energy  and  thermal  energy  correction  for  the  species  of  interest  are 
listed  in  Table  2.  The  calculated  gas  phase  AHf  (298K)  are  listed  in  Table  3  along  with  the  experi- 
mentally  measured  values  '  where  available  and  the  previous  calculated  BAC-MP4.  ’  ■  Table  4 
lists  the  gas  phase  AS0  (298K)  and  heat  capacity  (Cv)  calculated  by  using  the  HF/6-31g(d)  method. 


Table  2.  The  Electronic  Energy  Calculated  by  Using  MP2/6- 

31G(3df,2p)  Method  and  the  Thermal  Energy  Calculated  by 
Using  HF/6-31G(d)  Method  for  Species  of  Interest 


Species 

Electronic  energy  Hartree 

Thermal  correction  Hartree 

Si 

-288.88776 

0.00211 

Cl 

-459.6001 1 

0.0021 1 

H 

-0.4982329 

0.0021 1 

0 

-74.919245 

0.0021 1 

Ti 

-848.36909 

0.0021 1 

SiO 

-364.12262 

0.00582 

CM 

O 

0) 

-439.2289 

0.00993 

SiCI 

-748.6480 

0.00421 

SiClg 

-1208.4177 

0.00683 

SiClg 

-1668.1328 

0.01014 

SiCI4 

-2127.0469 

0.01345 

SiCIgOg 

-1818.2541 

0.01834 

SiCIgOgH 

-1818.9157 

0.03072 

SiCI20 

-1283.5642 

0.01188 

SiCIgOH 

-1743.9386 

0.02539 

TiO 

-923.5473 

0.00529 

Ti02 

-998.6500 

0.00684 

TiCI 

-1308.1312 

0.00403 

TiCI2 

-1767.9217 

0.00638 

TiClg 

-2227.6850 

0.00940 

TiCI4 

-2687.4362 

0.01255 

TiCIgOH 

-2303.4473 

0.02357 

TiCI20 

-1843,1225 

0.01093 

TiCIO 

-1383.3259 

0.00677 

TiCIgO 

-2302.7355 

0.01305 

TiCl302 

-2377.8071 

0.01737 
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Species  Electronic  energy  Hartree  Thermal  correction  Hartree 

TiOClOH  -1459.1236 

TiCIgOOH  -2378.4393 

Ti(OH)4  -1151.4594 


Table  3.  The  Gas  Phase  AHf°(298K)  Calculated  by  Using  Eqs.  3-5  as  Described 


in  the  Text.  Also  listed  are  previously  measured  experimental  value 
where  available  and  previous  calculated  values  using  BAC-MP4  method. 


Species 

AHj0(298K)  kj 
mol"1  This  work 

AHf°(298K)  k]  mol"1 
Previously  measured25 

AHf°(298K)  kj  mol"1 
BAC-MP415,20,21 

SiO 

-125.25 

-100.42 

-111.7 

SiC>2 

-361.87 

-305.43 

-282.00 

SICI 

151.60 

198.32 

157.99 

SiCI2 

-171.73 

-168.62 

-151.25 

SiClg 

-349.25 

-390.37 

-317.98 

SiCI4 

-701.72 

-662.75 

-662.70 

Si(OH)20 

-870.16 

-891 .2 

-804.1 

SiOHO 

-358.44 

-355.6 

-312.5 

SICI302 

-582.97 

-525.93 

SiCIgOgH 

-766.94 

-710.86 

SiCI20 

-511.52 

-465.26 

SiCIgOH 

-872.43 

-836.80 

TiO 

45.65 

54.39 

CM 

O 

H 

-188.15 

-305.43 

TiCI 

169.06 

154.39  (189.5)28 

TiCI2 

-208.72 

-237.23  (-196.6)28 

TiClg 

-513.53 

-539.32  (-497.5)ZB 

tici4 

-786.13 

-763.16 

TiCIgO 

-692.44 

-545.55 

TiCIO 

-303.29 

-244.26 

TiCIgO 

-604,93 

TiCl302 

-749.87 

TiOClOH 

-804.18 

TiCIgOOH 

-856.99 

Ti(OH)4 

-1274.56 

TiCIgOH 

-920.32 

0.02247 

0.02965 

0.06400 
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Table  4.  The  Gas  Phase  AS0  (298K)  and  Heat  Capacity  (Cv)  Calculated  by  Using 
HF/6-31g(d)  Method 


Species 

AS°(298K)  J  K'1  mol’1 

CV(298K)  J  1C1  mol'1 

AS0  (298K)  J  K"1  mol"1 
previously  measured25 

SiO 

21 1 .07 

21.22 

211.58 

Si02 

224.42 

33.06 

228.98 

SiCI 

238.15 

25.41 

237.81 

SiCI2 

280.44 

42.32 

281.30 

SiCI3 

325.58 

61.26 

318.13 

SiCI4 

347.54 

80.37 

330.86 

SiCIgC^ 

370.78 

93.36 

SiCI302H 

371.60 

100.96 

SiCIgO 

304.22 

55.34 

SiCIgOH 

346.82 

86.44 

TiO 

227.53 

21.83 

233.46 

CM 

O 

H 

297.12 

37.43 

260.14 

TiCI 

247.87 

27.04 

249.21 

TiClg 

228.66 

45.77 

278.34 

TiClg 

350.96 

66.57 

316.88 

TiCI4 

369.68 

86.41 

354.84 

TiCIgO 

327.21 

60.62 

TiCIO 

248.63 

32.40 

TiCIgO 

372.02 

83.80 

TiCIgC^ 

391.75 

100.59 

TiOClOH 

320.15 

67.43 

TiCIgOOH 

384.77 

104.39 

Ti(OH)4 

378.21 

125.82 

TiCIgOH 

389.89 

98.78 

Since  the  silicon  species  have  been  previously  calculated  using  BAC-MP4  and  many  have  been 
experimentally  measured,  this  allow  for  a  comparison  for  the  current  method  between  both  previous 
experimental  measurements  and  a  more  rigorous  theoretical  method  with  empirical  corrections.  The 
average  deviation  of  the  BAC-MP4  method  AHf°  (298K)  for  the  silicon  species  tested  compared  to 
previous  experimental  measurements  was  found  to  ±36.9  kJ  mol-1.  The  greatest  absolute  deviation 
from  experiment  was  for  SiCl3  and  Si(0H)20  with  a  deviation  of  72.4  kJ  mol-1  and  87.1  kJ  mol-1. 
Removing  these  species  reduces  the  average  deviation  to  ±22.6  kJ  mol-1.  The  current  method  pro¬ 
duces  an  average  deviation  of  the  BAC-MP4  method  AHf°  (298K)  for  the  silicon  species  tested  com¬ 
pared  to  previous  experimental  of  ±29.4  kJ  mol-1.  The  greatest  absolute  deviation  was  for  Si02  (56.4 
kJ  mol-"1)  although  SiCl3  once  again  has  one  of  the  greater  absolute  deviations  (46.7  kJ  mol-1). 
Removing  the  species  with  the  largest  deviation  slightly  reduces  the  average  deviation  to  ±25.5  kJ 
mol-1.  The  BAC-MP4,  in  general,  does  better  at  calculating  the  energies,  except  for  the  two  standout 
species  (SiCl3  and  Si(OH)20).  The  current  method,  however,  shows  that  it  does  an  adequate  job  of 
reproducing  the  AHf  (298K)  for  the  selected  species.  In  general  ,the  current  method  produces  a 
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slightly  more  negative  heat  of  formation  than  the  BAC-MP4  method  for  the  silicon  species  (on  aver¬ 
age  -38  ±  14  kJ  mor'lower). 

As  stated  above,  the  SiC,l3  AHf°  (298K)  has  a  large  deviation  from  the  experimental  method.  Irikura 
et  al.29  obtained  a  AHf°(298K)  of -353  ±  12  kJ  mol”1  by  combining  the  experimentally  measured 
AHf°  (298K)  of  SiCl3+  (41 8  ±  7  kJ  mol'1)30  with  an  ab  initio  calculation  of  the  IPa  of  SiCl  3  (7.93  ± 

0.05  eV).  This  heat  of  formation  is  in  better  agreement  with  the  BAC-MP4  and  current  calculations. 
Recent  measurements  of  Hildenbrand  et  al.  gave  a  AHf°  (298K)  of  SiCl3  &  -35 1 .5  kJ  mor1.31  They 
gave  a  recommended  value  of  AHf0  (298K)  =  -318  ±  4  kJ  mol'1  based  on  a  review  of  experiment  and 
theory.  This  value  is  in  better  agreement  with  the  current  predictions  of  both  this  study  and  the  BAC- 
MP4  study.  An  experimental  measurement  of  the  IPa  for  SiCl3  would  validate  a  lower  AHf°(298K) 
than  what  is  currently  recommended  by  NIST  for  SiCl3. 

The  BAC-MP4  method  predicts  the  heat  of  formation  of  SiO  and  Si02  accurately.  However,  it  does 
not  agree  with  the  current  measurement  of  the  Si(0H)20  heat  of  formation.  It  was  also  43  kJ  mol'1 
too  low  for  the  recently  updated  heat  of  formation  of  SiO(OH).  It’s  possible  that  the  Si-0  empirical 
correction  may  need  adjusting.  The  heat  of  formation  of  H(OH)SiO  has  previously  been  estimated  by 
group  additivity  to  be  -512  kJ  mol'1.32  BAC  MP4  predicts  this  heat  of  formation  to  be  -467.3  kJ 
mol'1.  The  current  method  can  be  used  to  predict  this  species  and  gets  a  result  slightly  closer  to  the 
group  additivity  method  (-495.7  kJ  mol'1).  Also,  the  bond  dissociation  energy  of  DH  (HSiOO-H)  = 
452  kJ  mol'1  was  used  by  Hartman  et  al.  This  dissociation  energy  is  ~82  kJ  mol'1  lower  than  what 
was  predicted  by  either  BAC  MP4  or  the  current  calculations. 

The  AHf  (298K)  of  the  titanium  species  was  found  to  have  a  larger  deviation  from  the  measured  val¬ 
ues  than  the  silicon  species.  This  was  not  unexpected  given  the  difficulty  that  results  in  treating  tran¬ 
sition  metals  by  ab  initio  methods.  The  current  method  produces  an  average  deviation  of  the  MP2 
method  AHf  (298K)  for  the  titanium  species  tested  compared  to  previous  experimental  values  of  ±53 
kJ  mol'1.  The  greatest  absolute  deviation  was  for  TiChO  (-146.9  kJ  mol'1)  and  Ti02  (1 17.3  kJ 
mol'1).  Removing  these  two  species  from  the  set  lowers  the  deviation  to  ±26.6  kJ  mol'1.  The  later 
number  was  more  typical  of  the  errors  seen  in  the  silicon  species.  Hildenbrand28  has  more  recently 
measured  the  AHf°(298K)  of  TiCl,  TiCl2  and  TiCl3  since  the  last  NIST  review  (189.5  kJ  mol'1, 

-196.6  kJ  mol'1  and  -497.5  kJ  mol'1),  respectively.  Using  these  values  produces  less  deviation  ±23.2 
kJ  mol'1. 

Similar  problems  were  observed  in  calculating  the  titanium  species  that  were  observed  in  the  silicon 
species.  The  silicon  species  with  the  highest  deviation  was  Si02,  and  TiO?  has  similar  problems. 
From  this,  it  seems  clear  that  this  level  of  theory  can  have  large  errors  associated  with  calculations  on 
close-shelled  species  that  contain  oxygen  double  bonds.  Although  several  closed-shell  species  with  a 
single  double  bond,  TiO  and  SiO  were  predicted  with  good  agreement  to  experimental  calculation. 

The  dissociation  energy  ( g)Q  at)  of  TiO  was  previously  measured  to  be  662.7  ±  6.7  kJ  mol'1,33  which 

is  in  good  agreement  with  our  calculation  (680  kJ  mol'1).  However,  the  Ti02  was  measured 

to  be  1271.9  ±  1 1.7  kJ  mol'1  (1260  ±  12  kJ  mol'1),  which  is  over  100  kJ  mol'1  off  from  the  calculated 
value  (1 162  kJ  mol'1).33,34  None  of  the  unknown  Ti  species  calculated  had  a  double-bonded  oxygen 
species.  Thus,  the  likely  error  associated  with  them  is  better  than  the  average  error  of  ±52.8  kJ  mol”1. 
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MP4  was  used  to  calculate  the  heat  of  formation  of  TiO,  Ti02,  SiO,  and  Si02  by  using  MP4/6- 
31G(3df,2p).  The  heat  of  formation  for  SiO  and  Si02  was  improved  by  using  MP4  (-107.63  kJ  mol"1 
and  -310.94  kJ  mol"1).  The  Ti02  was  also  slightly  improved  (-398.67  kJ  mol"1);  however,  the  TiO 
heat  of  formation  was  slightly  worse  (-23.81  kJ  mol"1).  MP2  can  provide  poor  geometries  for  some 
transition  metal  systems.  The  geometry  of  Ti02  was  optimized  by  using  density  functional  theory 
B3LYP  (Becke’s  three-parameter  exchange  functional  with  the  Yang-Parr  fit  for  the  correlation  func¬ 
tional)  instead  of  MP2.  However,  the  energy  obtained  was  not  an  improvement  on  the  original 
method 

The  calculated  AS0  (298K)  for  the  silicon  species  was  found  to  have  an  absolute  deviate  from  experi¬ 
ment  by  only  ±5  J  K"1  mol"1.  The  calculation  for  the  Ti  species  was  found  to  be  worse  with  a  devia¬ 
tion  of  23  J  K"1  mol"1.  Overall,  the  method  does  a  good  job  at  predicting  the  entropy  of  the  species  of 
interest. 
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4.  Conclusions 


Titanium  intermediates  associated  with  the  combustion  of  TiCU  in  an  O2  flame  have  been  calculated 
by  using  MP2/6-31G(3df,2p).  The  theory  does  an  adequate  job  of  reproducing  known  heats  of  for¬ 
mation  for  silicon  and  titanium  compounds.  Except  for  titanium  compounds  that  contain  four  bonds 
total  with  two  of  them  being  an  oxygen  double  bond.  Empirical  corrections  will  be  needed  for  these 
Ti=0  species  in  order  to  adequately  predict  them  using  ab  initio  methods. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security  programs,  specializing 
in  advanced  military  space  systems.  The  Corporation's  Laboratory  Operations  supports  the  effective  and 
timely  development  and  operation  of  national  security  systems  through  scientific  research  and  the  application 
of  advanced  technology.  Vital  to  the  success  of  the  Corporation  is  the  technical  staffs  wide-ranging  expertise 
and  its  ability"  to  stay  abreast  of  new  technological  developments  and  program  support  issues  associated  with 
rapidly  evolving  space  systems.  Contributing  capabilities  are  provided  by  these  individual  organizations: 

Electronics  and  Photonics  Laboratory:  Microelectronics,  VLSI  reliability,  failure  analysis, 
solid-state  device  physics,  compound  semiconductors,  radiation  effects,  infrared  and  CCD 
detector  devices,  data  storage  and  display  technologies;  lasers  and  electro-optics,  solid-state 
laser  design,  micro-optics,  optical  communications,  and  fiber-optic  sensors;  atomic  frequency 
standards,  applied  laser  spectroscopy,  laser  chemistry,  atmospheric  propagation  and  beam 
control,  LIDAR/LADAR  remote  sensing;  solar  cell  and  array  testing  and  evaluation,  battery 
electrochemistry,  battery  testing  and  evaluation. 

Space  Materials  Laboratory:  Evaluation  and  characterizations  of  new  materials  and 
processing  techniques:  metals,  alloys,  ceramics,  polymers,  thin  films,  and  composites; 
development  of  advanced  deposition  processes;  nondestructive  evaluation,  component  failure 
analysis  and  reliability;  structural  mechanics,  fracture  mechanics,  and  stress  corrosion; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch  vehicle 
fluid  mechanics,  heat  transfer  and  flight  dynamics;  aerothermodynamics;  chemical  and  electric 
propulsion;  environmental  chemistry;  combustion  processes;  space  environment  effects  on 
materials,  hardening  and  vulnerability  assessment;  contamination,  thermal  and  structural 
control;  lubrication  and  surface  phenomena  Microelectromechanical  systems  (MEMS)  for 
space  applications;  laser  micromachining;  laser-surface  physical  and  chemical  interactions; 
micropropulsion;  micro-  and  nanosatellite  mission  analysis;  intelligent  microinstruments  for 
monitoring  space  and  launch  system  environments. 

Space  Science  Applications  Laboratory:  Magnetospheric,  auroral  and  cosmic-ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using  atmospheric 
radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis;  infrared  surveillance, 
imaging  and  remote  sensing;  multispectral  and  hyperspectral  sensor  development;  data 
analysis  and  algorithm  development;  applications  of  multispectral  and  hyperspectral  imagery 
to  defense,  civil  space,  commercial,  and  environmental  missions;  effects  of  solar  activity, 
magnetic  storms  and  nuclear  explosions  on  the  Earth’s  atmosphere,  ionosphere  and 
magnetosphere;  effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation,  design,  fabrication  and  test;  environmental  chemistry,  trace  detection; 
atmospheric  chemical  reactions,  atmospheric  optics,  light  scattering,  state-specific  chemical 
reactions,  and  radiative  signatures  of  missile  plumes. 


